Optical interferometry in the presence of large phase diffusion 
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Phase diffusion represents a crucial obstacle towards the implementation of high precision interferometric 
measurements and phase shift based communication channels. Here we present a nearly optimal interferometric 
scheme based on homodyne detection and coherent signals for the detection of a phase shift in the presence of 
large phase diffusion. In our scheme the ultimate bound to interferometric sensitivity is achieved already for a 
small number of measurements, of the order of hundreds, without using nonclassical light. 
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I. INTRODUCTION 



Optical interferometry represents a high accurate measure- 
ment scheme with wide applications in many fields of sci- 
ence and technology fV-'5]. Besides, the precise estimation 
of an optical phase shift is relevant for optical communica- 
tion schemes where information is encoded in the phase of 
travelling pulses. Several experimental protocols have been 
proposed and demonstrated to estimate the value of the op- 
tical phase lIB Um and showing the possibility to attain the 
so-called Heisenberg limit 1, 1214221 . Recent developments 
also revealed the potential advantages of nonlinear interac- 
tions |23|. However, in realistic conditions, one has to re- 
trieve phase information that has been unavoidably degraded 
by different sources of noise, which have to be taken into ac- 
count in order to evaluate the interferometric precision [24]. 
The effects of imperfect photodetection in the measurement 
stage, or the presence of amplitude noise in the interferomet- 
ric arms have been extensively studied 1251 435 1 . Only recently, 
the role of phase-diffusive noise in interferometry have been 
theoretically investigated for optical polarization qubit ll36l - 
|38J . condensate systems 1,39. ,40] . Bose-Josephson junctions 
ll4n . and Gaussian states of Ught B2ll . As a matter of fact, 
phase-diffusive noise is the most detrimental for interferome- 
try and any signal that is unaffected by phase-diffusion, is also 
invariant under a phase shift, and thus totally useless for phase 
estimation. 

In this paper, we present an experimental interferometric 
scheme where phase diffusion may be inserted in a controlled 
way, and demonstrate that homodyne detection and coher- 
ent signals are nearly optimal for the detection of a phase 
shift in the presence of large phase diffusion. Indeed, while 
in ideal conditions squeezed vacuum is the most sensitive 
Gaussian probe state for a given average photon number ll43l . 
for large phase-diffusive noise, coherent states become the 



optimal choice, outperforming squeezed states (42]. In our 
scheme the ultimate bound to interferometric sensitivity, as 
dictated by the Cramer-Rao (CR) theorem, is achieved already 
for a small number of repeated measurements, of the order of 
hundreds, using Bayesian inference on homodyne data and 
without the need of nonclassical light. 

The paper is structured as follows: In Section [ll] we de- 
scribe the evolution of a light beam in a phase diffusing envi- 
ronment as well as the bound to interferometric precision in 
the presence of phase noise. In Section [III] we describe our 
experimental apparatus, whereas the experimental results are 
reported and discussed in Section IV SectionMcloses the pa- 
per with some concluding remarks. 



n. INTERFEROMETRY IN THE PRESENCE OF PHASE 
DIFFUSION 

The evolution of a light beam in a phase diffusing envkon- 
ment is desciibed by the master equation 







rCla'alg, 



where C[0]g = 20gO^ - O^Og - gO^O and T is the phase 
damping rate. An initial state go evolves as 



ft =AAA(eo) ^E^"^'^"""^'^".™! 



n){m\ , 



where A = Ft, and £»„_,„ = {n\go\m). The diagonal ele- 
ments are left unchanged, in fact energy is conserved, whereas 
the off-diagonal ones are progressively destroyed, together 
with the phase information carried by the state. Phase dif- 
fusion corresponds to the application of a random, zero-mean 
Gaussian-distributed phase shift, i.e., 

gt = j^d(3g{(3\A)UpgoUl g(/3|A) = — ^=^ (1) 



•Electronic address: |matteo.paris@fisica.ununi.it| 



where Up = exp{— i/3(a' a)} is the phase shift operator. 



We assume that the phase noise occurs between the apph- 
cation of the phase shift and the detection of the signal, and 
consider the estimation of a phase shift appHed to a single- 
mode coherent state. Homodyne detection is then performed 
on the output state 

£'A.a((/>) = NA{U^\a){a\Ul) , 

and the value of the unknown phase shift </) is inferred us- 
ing Bayesian estimation applied to homodyne data. Notice, 
however, that since the phase noise map and the phase shift 
operation commute, our results are valid also when the phase 
shift is applied to an already phase-diffused coherent state. 
The precision of the above procedure is then compared with 
the benchmarks given by i) the quantum CR bound for co- 
herent states and any quantum limited kind of measurement, 
ii) the ultimate precision achievable with optimized Gaussian 
states, i.e., the quantum CR bound for general Gaussian sig- 
nals, where, e.g., we allow for squeezing. 

A. Interferometric precision in the presence of phase noise 



The quantum CR bound i PTTIPi SI is obtained starting from 
the Born rule p(x|0) = Trpa;^!^] where {Il^r} is the operator- 
valued measure describing the measurement and g^ the den- 
sity operator of the family of phase-shifted states under in- 
vestigation. Upon introducing the (symmetric) logarithmic 
derivative L^ as the operator satisfying 



where 
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Assuming that no a priori information is available on the value 
of the phase shift (i.e., uniform prior), and using the Bayes 
theorem, one can write the a posteriori probability 



P(cj,\X)^^L{X\cj,) 



N- 



d(t>L{(t>\X), (2) 



$ = [0, tt] being the parameter space. The probability 
P{(I)\X) is the expected distribution of (j) given the data sample 
X. The Bayesian estimator 0b is the mean of the a posteriori 
distribution, whereas the sensitivity of the overall procedure 
corresponds to its variance 



= / d(A(0-</)B)'P(<^|X) 
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Bayesian estimators are known to be asymptotically unbiased 
and optimal, namely, they allow one to achieve the CR bound 
as the size of the data sample increases B9ll50l . On the other 
hand, the number of data needed to achieve the asymptotic 
region may depend on the specific implementation 1 5 1 1 . In the 
following we will experimentally show that our setup achieves 
optimal estimation already after collecting few hundreds of 
measurements. 



one proves that the ultimate limit to precision (independently 
on the measurement used) is given by the quantum CR bound 

Var(0) > [MH{(t))]-\ 

where H{(j)) ~ Tr[p^ LV\ is the quantum Fisher information 
(QFI). The ultimate sensitivity of an interferometer thus de- 
pends on the family of signals used to probe the phase shift 
and thus, as said above, we are going to compare the precision 
of our interferometer with the maximum achievable with co- 
herent states, and with the ultimate precision achievable with 
optimized Gaussian states (for more details about the deriva- 
tion of the corresponding quantum CR bounds see ll42l ). 
Homodyne detection measures the field quadrature 
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where 9 — arg a + 7r/2 is set to the optimal value to detect 
the imposed phase shift. The likelihood of a set of homodyne 
data 

X = {xi,X2,.--,Xm} , 

is the overall probability of the sample given the unknown 
phase 4>, i.e., 

M 



LiX\cj,) = l[p{xk\q^), 
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ni. EXPERIMENTAL APPARATUS 

A schematic diagram of the interferometer is reported in 
Fig.[T] The principal radiation source is provided by a He:Ne 
laser (12 mW, 633 nm) shot-noise limited above 2 MHz. The 
laser emits a linearly polarized beam in a TEMOO mode. The 
beam is splitted into two parts of variable relative intensity 
by a combination of a halfwave plate (HWP) and a polariz- 
ing beam splitter (PBS). The strongest part is sent directly to 
the homodyne detector where it acts as the local oscillator, 
whereas the ramaining part is used to encode the signal and 
will undergo the homodyne detection. The optical paths trav- 
elled by the local oscillator and the signal beams are carefully 
adjusted to obtain a visibility typically above 90% measured 
at one of the homodyne output ports. The signal is amplitude 
modulated at 4 MHz with a defined modulation depth to con- 
trol the average number of photons in the generate state. 

The amplitude modulation system consist of a KDP non- 
linear crystal with the xy axes at 45°, and a PBS. The mod- 
ulation is applied at the KDP crystal by means a waveform 
generator Rohde & Schwarz and a power amplifier Mini- 
Circuits ZHL-32A. The modulation depth is imposed at the 
proper level by a computer that sends a costant voltage to a 
mixer (Ml) located between the waveform generator and the 
power amplifier. One of the mirrors in the signal path is piezo 
mounted to obtain a variable phase difference between the two 
beams. The piezo is preloaded and its resonance frequency is 
13.5 kHz. 



The phase difference is controlled by the computer after a 
calibration stage. The computer sends a voltage signal be- 
tween and 10 V that corresponds at the phase diffusion with 
a frequency of 5 kHz to a power amplifier based on LM675 in- 
tegrated circuit that is able to drive the piezo at this frequency. 
With this system it is possibile to generate any kind of phase 
modulation. 



HWP PBS 




FIG. 1: (color online). Schematic diagram of the experimental setup. 
A He:Ne laser is divided into two beams, one acts as the local oscil- 
lator and the other represents the signal beam. The signal is modu- 
lated at 4MHz with a defined modulation depth to control the aver- 
age number of photons in the generate state. One of the mirrors in 
the signal path is piezo mounted to obtain a variable phase difference 
between the two beams. The data are recorded by a homodyne detec- 
tor whose difference photocorrent is demodulated and then acquired 
by a computer after a low pass filter. We also show the typical homo- 
dyne samples obtained for coherent signals of different amplitudes 
by varying the phase of the local oscillator (these are used to check 
the calibration of the piezo, which is performed using signals with a 
larger number of photons). 

The detector is composed by a 50:50 beams splitter (BS) 
and a balanced amplifier detector with a bandwidth of 
50 MHz. The difference photocurrent is filtered with high pass 
filters, amplified and demodulated at 4 MHz by means of an 
electrical mixer (M2). In this way the detection occurs outside 
any technical noise and, more importantly, in a spectral region 
where the laser does not carry excess noise. The signal is fil- 
terd by a low pass filter with a bandwidth of 300 kHz and sent 
to the computer through the National Instrument multichannel 
data acquisition 6251 with 16 bit of resolution and 1.25 MS/s 
sampling rate. The same device is used to send diffusion pa- 
rameters to the phase modulator and signal parameters to the 
amplitude modulator. 



IV. EXPERIMENTAL RESULTS 

In this Section, we present our experimental results, ob- 
tained with signals of different energies and different levels 
of noise. At first we show homodyne samples with the corre- 
sponding a posteriori distributions and then compare the pre- 
cision obtained in our scheme with the ultimate bound im- 
posed by the (quantum) Cramer-Rao theorem. Finally, we 
analyze the dependence of precision on the signal energy 



and the noise in order to illustrate how in the limit of large 
phase diffusion coherent states becomes the optimal Gaus- 
sian probe states. In fact, they outperform squeezed vacuum 
states, whose non-classical features are degraded by phase dif- 
fusion process, to an extent that make them useless for quan- 
tum metrology. 

In Fig.l2]we report typical examples of homodyne samples, 
referred to a coherent signal with N = |ap mean photon 
number measured at fixed optimal 9, together with the cor- 
responding Bayesian a posteriori distribution for the phase 
shift. The yellow area denotes the portion of data used to 
infer the phase shift. We choose this range in order to em- 
phasize that the optimality region in achieved already in that 
region. In fact, upon considering larger samples, precision 
would be improved, due to the statistical scaling of the vari- 
ance Var[0] = C/M, C being a proportionality constant. On 
the other hand, optimality, i.e., the fact that 

C - 1/H^ , 

where Ha is the QFI for phase-diffused coherent signals, is 
achieved for M ^ 100 measurements. In the noiseless case 
the QFI is given by Ha — 4:N, whereas it decreases monoton- 
ically by increasing the value of the noise parameter A. No- 
tice that using optimized Gaussian signals, i.e. the squeezed 
vacuum state, one has a QFI given by Hg — 8iV^ + 8N in 
the noiseless case. However, in the presence of large phase 
diffusion, i.e. for large values of A, Ha is larger than the QFI 
obtained for phase-diffused squeezed vacuum states. In other 
words, coherent states turns out to be the optimal Gaussian 
probe states II42I . 
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FIG. 2: (color online) Typical examples of homodyne samples mea- 
sured at fixed optimal 6, together with the corresponding Bayesian 
a posteriori distribution for the phase shift. The phase diffusion is 
A — 7r/6 rad and the yellow area denotes the portion of data used to 
infer the phase shift. 



In Fig.|3]we plot the quantity 

i^M = MVar[0B]i?„, 

i.e., the variance of the Bayesian estimator from homodyne 
data multipHed by the number of data (measurements) and by 
the coherent states quantum Fisher information, as a function 
of M. Km is by definition larger than one and expresses the 
ratio between the actual precision of the interferometric setup 
and the CR bound. As it is apparent from the plot Km rapidly 
decreases with the number of measurements, almost indepen- 
dently on the value of the number of photons N and of the 
noise parameter A. The optimality region, i.e.. Km — 1 is 
achieved already for M ~ 100 measurements, and the asymp- 
totic value of Km is closer to 1 for increasing N and A. Fur- 
thermore, the number of measurements needed to achieve the 
optimal region may be (slightly) reduced by using the Jeffreys 
prior 1,5 2 J 



instead of the uniform one, where 



F((j))^ j dxp{x\(j))[d^\ogp{x\ 



is the Fisher information of the homodyne distribution. 
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FIG. 3: (color online) The noise ratio Km ~ (Va,r[(j}B]MHa) as 
a function of the number of data AI and for different values of the 
number of photons A'^ and the noise parameter A. Blue circles: TV = 
0.90, A = 7r/18 rad; red squares: iV = 0.90, A = 7r/9 rad; yellow 
diamonds: A'' = 4.12, A = 7r/18 rad; green triangles: A^ = 4.12, 
A = 7r/9 rad. 

In Fig. [4] we show the variance of the Bayesian estimator 
from homodyne data 

Vm = MVar[0B] 



obtained after M measurements, together with the CR bound 
1/Ha for coherent states, and for the (phase-diffused) opti- 
mized Gaussian states, i.e., l/Hg. In particular, the top panel 
shows the behaviour as a function of A for different values of 
the number of photons N, while in the bottom panel we plot 
the same quantities as a function of the number of photons N 
and for different values of the noise A. As it is apparent from 
the plots, nearly optimal inferferometric precision is achieved 
for increasing energy or phase diffusion, i.e., for larger values 
of N or A. 



V. CONCLUSIONS 

In conclusion, we have demonstrated a nearly optimal in- 
terferometric scheme based on homodyne detection and co- 
herent signals for the detection of a phase shift in the pres- 
ence of large phase diffusion. Our scheme does not require 
nonclassical light and achieve the ultimate bound to interfer- 
ometric sensitivity using Bayesian analysis on small samples 
of homodyne data, where the number of measurements is of 
the order of few hundreds. 

It is worth noting that for large phase diffusion coherent 
states are the optimal Gaussian probe states. Indeed they out- 
perform squeezed vacuum states, whose non-classical features 
are degraded by phase diffusion process, such that they be- 
come completely useless for quantum metrology. 

Optical interferometry represents a high accurate measure- 
ment scheme with wide applications in many fields of science 
and technology, including high precision measurements and 
communication channels. On the other hand, phase diffusion 
represents a crucial obstacle towards the implementation of 
high precision interferometric measurements and phase shift 
based communication channels. Our results allow to design 
feasible, high-performance, communication channels also in 
the presence of phase noise, which cannot be effectively con- 
trolled in realistic conditions. Therefore, besides fundamental 
interest, our results also represent a benchmark for realistic 
phase based communication or measurement protocols. 
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